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1. Introduction 

The asymptotic expansion of the heat kernel proves to be important in a series of 
physical applications. Its coefficients specify the divergences and conformal anomalies 
taking place in a concrete field theory model the high temperature behaviour of the 
thermodynamic functions P, and so on 0]. 

For a flat manifold with a smooth boundary the heat kernel coefficients are 
determined by local characteristics of the boundary |5], || [7|, ||. If the boundary has 
discontinuities (for example, it is piecewise smooth), then the latter give additional 
contributions to the heat kernel coefficients [f|. Usually in physical applications one 
assumes that the boundary is smooth. However there is a series of problems where such 
assumption is certainly not acceptable. A typical example here is supplied by fields 
defined between two plates which cross at a given angle, i.e., inside a dihedral angle. 
Such a configuration is considered when calculating the Casimir effect for a conducting 
wedge |1J. If the fields inside a dihedral are subjected to periodicity condition with 



respect to the angular variable then one is concerned in fact with the fields on a cone, 
and the point where the boundary has discontinuity becomes an internal point of the 
cone surface. However the origin of this singularity is the same as in the case of fields 
inside the wedge. The conical singularity proved to be very important in many areas of 
mathematical physics |TI, [12] and lately it is investigated in connection with studies of 



quantum fields on the background of black holes |TB[] and cosmic strings [[H 



The general consideration of the boundary nonsmothness in terms of the heat kernel 
expansion lacks till now \TE\. Such contributions to the heat kernel coefficients C3/2 



and C2 have been investigated in papers M. For a plane domain the contribution 



to C\ generated by an edge of the boundary is known [|Ty, [17| and by its limiting 
configuration, by cusp. It is interesting to note that a cusp pointing outward, with 
respect to the domain under study, leads to change of the power of the asymptotic 
variable (time) in the heat kernel expansion || in comparison with the standard case. 
In references [||, [Tj| [18|, [19| the asymptotic expansion of the heat kernel with allowance 



for the boundary nonsmothness has been built by calculating the relevant Green function 
of the heat conduction equation. The present paper seeks to show the effectiveness of 
applying the spectral zeta functions for the calculation of the contributions to the heat 
kernel coefficients caused by such boundary discontinuities as the corners. For this goal 
we shall use the technique for constructing the zeta functions developed in [[2(J and 



extended in |[21||. This method proves to be very effective for calculating the heat kernel 



coefficients for different boundary conditions given on a sphere and cylinder || 20, 22] 



A close approach has been used in E3L 24 



We shall consider internal and external parts of a plane sector formed by two infinite 
radial rays emerging from the center of a circle of radius R at angle a to one another 
(see figure 1, where I is the internal circular sector and II is the external circular sector). 
The choice of such domains with nonsmooth boundary is caused by the possibility of 
constructing for them the global zeta functions. The latter cannot be done, for example, 
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for an open angle (figure 1 without circular arc inside the angle). This point will be 
discussed in detail in Section II. In both the sectors the Laplace operator is defined acting 
on scalar real functions subjected to the Dirichlet or Neumann boundary conditions. If 
we substitute in figure 1 the radial rays by crossed planes and the circular arc 1-2 
by an appropriate part of a cylinder surface we arrive at the boundary value problems 
which have the same heat kernel coefficients as in the plane case. 

Identifying the points of the boundary with the polar coordinates (r, 9 = 0) and 
(r, 9 = a), i.e., imposing the periodicity condition with respect to the angular variable 
9 with a period a, we arrive obviously at the spectral problem for the Laplace operator 
on two parts of lateral surface of a cone C a : internal part (r < R) and external part 
(r > R). At r = R we can, as before, impose the Dirichlet or Neumann conditions. 

In the present paper, six coefficients of the asymptotic expansion of the heat kernel 
for the boundary value problems specified above will be calculated by making use of the 
relevant zeta functions. It will be shown that each of these coefficients, starting from 
the third one, is the sum of contributions generated by the corners of the boundary 
and by the curvature of the arc 1-2. Analysis of the obtained results enables one to 
reveal some regularities obeyed by the contributions of the boundary nonsmoothness to 
the heat kernel coefficients, namely: i) the corner contribution substantially depends on 
whether the corner is made up by crossing two straight lines or two lines with nonzero 
curvature; ii) contributions of adjacent angles to the heat kernel coefficients B 3 / 2 and 
B 5 / 2 have opposite signs in the same way as the contributions of a circular arc to the 
coefficients B\ and B 2 for the internal and external regions. In the case of a polygon or 
its cylindrical generalization the rules are formulated for obtaining all the heat kernel 
coefficients for the minus Laplace operator. The implications of the obtained results 
in the Casimir energy calculations with employment of the zeta function technique are 
also considered. In particular, it is shown that nonsmoothness of the boundary does not 
make always worse the situation with calculation of the vacuum energy in the framework 
of the zeta regularization method. For example, the corner contributions can simply 
be mutually cancelled in the same way as the contributions of the curved boundary 
are cancelled when taking into account the internal and external regions. Besides, 
the analysis conducted enables one to infer that it is very unlikely to get a finite and 
unique result for the vacuum energy by smoothing the boundary singularities, e.g., by 
taking into account the atomic structure of the boundary or quantum fluctuations of 
the boundary. 

The paper is organized as follows. Section 2 is devoted to the detailed discussion 
of the choice of the domains with piece-wise smooth boundary (internal and external 
circular sectors), for which the complete spectral zeta functions can be constructed. In 
section 3 the heat kernel coefficients are calculated for internal circular sector by making 
use of the counter integral representation for the corresponding zeta functions ^(| |[. 
In section 4 the contributions to the heat kernel coefficients of the individual boundary 
discontinuities (corners of different angles) are identified. Further (in section 5) the heat 
kernel coefficients for the union of both the sectors are calculated and the coefficients 
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for the external sector are obtained as the regarding differences. Technically it turns 
out to be simpler in comparison with calculation of the heat kernel coefficients for the 
external sector alone. The asymptotic expansion of the heat kernel for the union of 
internal and external sectors is constructed by differentiation of the logarithm of the 
function specifying the frequency equation [f23|, |24]]. In section 6 the identification of 
the individual contributions of boundary nonsmoothness to heat kernel coefficients for 
external sector is carried out. Some general rules for these contributions are revealed 
here for both sectors. In section 7 we conclude with a few summarizing remarks. 



2. Choice of domains with piecewise smooth boundaries 

When choosing the domain for investigating the contribution to the heat kernel 
coefficients of the boundary discontinuities we pursue two goals: the boundary of a 
domain should have a sufficient number of discontinuities and at the same time for this 
domain one can construct the spectral zeta function. In order to investigate the corner 
singularities of the boundary one should at first sight take the most simple configuration, 
namely, the angle on plane formed by two radial unrestricted rays or dihedral in space. 
However, for these domains the global zeta functions cannot be constructed. Let us 
explain this point in detail. We consider the Laplace operator A acting on the scalar 
functions defined inside the dihedral of opening angle a (wedge of an angle a, W a x M 1 ) 
and subjected to the Dirichlet conditions on the wedge sides. In cylindrical coordinate 
system (r, 9, z) the eigenfunctions in this problem are 

u X nk(r,9,z) = —=J np (\r)sm(np9), 



0<9<a, p = n/a, n=l,2,..., < A < oo , (2.1) 

where J v {z) is the Bessel function. The operator —A has the following eigenvalues 

u 2 {k, A) = k 2 + A 2 , -oo < /c < oo, < A < oo . (2.2) 

These eigenvalues do not depend on the quantum number n, i.e., there is a degeneracy 
with respect to this number, the multiplicity of this degeneracy being infinite^ 

oo 

iV = ^l = oo. (2.3) 

71=1 

The global spectral zeta function of the operator L, C( s ), is determined as the trace 
of the operator L~ s 

COO = ^ L-* = J^-A;* (2.4) 

j 

§ In addition to infinite multiplicity, the every point of the spectrum with fixed values of k and A is a 
nonisolated point. In fact, the spectrum is continuous according to the terminology of spectral theory 
of operators in Hilbert space |^5|, and tu 2 (k, A) are not eigenvalues as the eigenfunctions are not 

square integrable. 
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where A| is the jth eigenvalue of the operator L, and Nj is the degeneracy of this 
eigenvalue. Obviously, this definition does not work when Nj is infinite. 

In order to remove this degeneracy we put inside a dihedral a cylindrical boundary 
as it is shown in figure 1 (the arc 1-2). On the internal and external sides of this 
boundary the scalar field will obey the Dirichlet or Neumann conditions. Thus, we are 
considering the internal (I) and external (II) sectors. Certainly, the union of these two 
sectors is not identical to unrestricted dihedral (to a wedge), because now the values of 
the field on the arc 1-2 are not arbitrary but they are determined by the corresponding 
boundary conditions. In the new configuration there appear additional discontinuities 
of the boundary at the points 1 and 2. However at these points the angle, at which 
the involved boundary surfaces intersect, is fixed (it is equal to vr/2). It will be shown 
below that the contribution of such boundary singularities to the heat kernel coefficients 
can be easily separated from the contribution of the corner at the origin. The latter 
is proportional to the difference n — a, because at a = ir the singularity at the origin 
disappears. 

We shall concern with the standard asymptotic expansion of heat kernel 

oo 

K(t) = e ~ X2]t = i^t)- d/2 tnl2R n/2 + ES, (2.5) 

j n=0 

where d is the dimension of the manifold under study and ES stands for the exponentially 
small corrections as t — > 0. This definition leads to the same heat kernel coefficients 
B n /2 for a dihedral and for corresponding plane problem obtained by crossing of the 
dihedral by a transverse plane. In fact, the eigenvalues of the operator (—A) in these 
two problems obey the relation 

\i(d = 3) = k 2 + \ 2 J (d = 2). (2.6) 

Hence 

K d=3 (t) = ^e- kH K d=2 {t) = —L= K d=2 {t). (2.7) 

Taking into account the definition fl2.5|) one easily deduce from fl2.7|) that the heat kernel 
coefficients B n / 2 in two eigenvalue boundary problems, mentioned above, are equal. Here 
one should bear in mind that the heat kernel K d=3 (t) and its coefficients are referred 
to a unite length along the OZ axes. In view of this, when calculating the heat kernel 
coefficients we shall consider either the spectral problem on a plane (d = 2) or in the 
space (d — 3), pursuing only simplicity of calculation. 

3. Heat kernel coefficients for internal sector 

3.1. Dirichlet boundary conditions 

At first we consider only internal sector (the region I in figure 1) for d = 3, i.e., for 
internal sector of a dihedral or a wedged]]. We employ here the technique developed 

| In the framework of the quantum billiard studies this problem has been investigated in fl^| 
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in |Mj and extended to the generalized bounded cone in |2T[ (see also book ||). In this 
approach the spectral zeta function ((s) should be constructed in the beginning, and 
then the relevant heat kernel coefficients are calculated through the relation || 

lim (s + n-d/2)((s)r(s), n = 0, 1/2, 1, . . . . (3.1) 



(47l) d / 2 s^d/2-n 

For the Dirichlet boundary conditions the eigenfunctions of the operator (—A) in 
the region I are defined by equation ( j2.1p , with A being the roots of the equations 

J np {KmR) = 0, p = n/a, n = l, 2, .... (3.2) 

Here the subscript m = 1, 2, . . . numbers the nonzero roots of these equations for fixed 
n. The relevant eigenvalues are defined by equation ( |2.2|) with A = X nm 

u 2 = k 2 + A^„, — oo < k < oo . (3.3) 

In view of the behavior of the Bessel function near zero J v (z) ~ z u /(2 u T(y+ 1)) it follows 
that the frequency equation (|3.2| ) has the zero root of 'multiplicity' v. Such roots should 
be removed from the definition of the spectral zeta function (|2.4j) . Therefore instead of 
equation (|3~2^ ) we shall use the following frequency equation 

(XR)- U J V (XR) = 0, u = np, n = l,2, ... . (3.4) 

According to the general definition (2.4) the zeta function in the problem under study 
is given by 

^EE( fc2 + A nJ- S - (3.5) 

As usual we substitute the sum over m in this formula by the contour integral in the 
plane of a complex variable A 

Cd{s) = r § £ -hi f dA(A2 + k2) ~ s l\ i ^ XR y uj ^ XR ^ ' ( 3 - 6 ) 

J -OO n _y JC 

where the contour C encloses counter clockwise the positive roots of equations ( |3.4|) . 
On deforming the contour C in an appropriate way and integrating with respect to k 
by means of the formula 

d " -V^ r( *-y 2) , Re S >i, (3.7) 



C D (s) = } R / P) .„ ,f n 1 - 28 / °dyy 1 - 2s ^ln[(z/ 2 /)-%(^)],(3.8) 



(P + A 2 ) s v r(s) ' 2' 

we arrive at the result 

) 2s 

2 vfirooFg - a ) ^ 

where /^(^) is the modified Bessel function. 

The analytical continuation of ( p.8[ ) into the left half-plane of the complex variable 
s is accomplished by making use of the uniform asymptotic expansion of the function 

Iviyy) 
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m — 2,1 y dr ? v 1 + y (o-,n\ 
v = v 1 + y + ln ~ n — = . -r = - > 3 - 10 

1 + y/l + y 2 dy y 



and Uk(t) are the known polynomials in t = + y 2 - Their explicit form and the 

corresponding recurrent relations can be found in references [28, |29[| . Keeping in this 
expansion all the terms proportional to v~ k with k < 4 one can write 

where 





k=l J 








1 5 , 

= -t 1 3 , 

8 24 ' 
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16 ' 
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«= 221 7 
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In accordance with this expansion we represent the zeta function in equation (|3.9| ) 
as the following sum 

4 

where 

/■OO J 

Z? 1 (s) = C(s)p l - 2 'C R (2s-2) dyy 1 " 2 *— fo-lny), (3.14) 

Jo dy 

= - l C(s)p 1 - 2 Xn(2s~l) J dyy^-Hl+y 2 ), (3.15) 

f°° H 

Zf(a) =C( S y- 2s - 1 U2s + j-l) dyy^-F^t), j = 1,2,3,4. (3.16) 
In these formulas Cr( s ) is the Riemann zeta function and 

CW = v?r(.W-,) - (3 ' 17) 

Substituting (|3.10|) into equation Q3.14Q and taking into account the value of the 
integral 

d „ v 2 - 2 ' _i r( 3 -i)r(3/2- S ) 

we obtain for the function Z^^s) the final expression 

Z° l ( s ) = -^=(-\ 2S 2 - r( '~ 1) / , ( R (2s-2). (3.19) 
~ lW 8v^VP/ (s-1/2)1» S/U ; 1 ; 
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The product T(s) Z^s) has simple poles at the points s = 3/2, 1, and 1/2 with the 
respective residua 

aR 2 R 



16vr 3 / 2 ' 8tt' 48 a' (3 ' 20) 
Multiplying these residua by (47r) 3 / 2 we obtain, according to equation (|3.1|) , the 
contributions of the function Z^^s) to the heat kernel coefficients B , -B1/2, and £>i, 
respectively. 



By making use of the table integral [ 30 1 



l + 



2/ T^ = r(^- s )r( s -±), (3.21) 



^-^S ^-l)^. (3-22) 



we recast equation ( |3.15| ) to the form 

2s-l 

'~Z \p ) ' ' I'j.s 

The expression T(s) (s) has simple poles at the points s — 1/2 and s = 1 with the 
residua 



16 tt' 16^' (3 ' 23) 

which give the respective contributions to the heat kernel coefficients Bi/ 2 and Bi (see 
table 1). 

Calculation of the functions Zj(s), j — 1, . . . , 4, defined in equation (|3.16|) , can be 



carried out with hardly any trouble. Their contribution to the heat kernel coefficients 
are given in table 1. In order to get the complete values of these coefficients one should 
sum all the elements of respective rows. 

3.2. Neumann boundary conditions 

In this case the operator —A has the following eigenfunctions in internal circular sector 
I (see figure 1) 

e ikz 

v\nk(r,0,z) = r] n0 —= J np (X nm r) cos(np6), n = 0, 1,2, . . . , 
\/ira 



(3.24) 



Here X nm are the roots of the equations 

J' np (X nm R) = 0, n = 0, 1, 2, . . . . (3.25) 

Taking into account the behavior of the derivative of the Bessel function at the origin 
J' u (z) ~ z u - 1 /(2 u T(u)), we multiply equation ( ^2"5l) by (\ nm R) 1 ~ np in order to exclude 
the zero multiple roots 

(A Rf- V 'J' V (XR) = 0, v = np, n = 0,l, ... . (3.26) 
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By making use of the frequency equations ( |3.26| ) we can write immediately the 
integral representations for the zeta function in the problem at hand analogous to 
equation (|3.8|) 



Cn(s) = C(s) 



l-2s 



d_ 

dy 



Hyi'oiy)) 



(3.27) 



°° POO J 

+ 5> 1 - 2s / dyy^-^-Hy^lUuy)) 

n=l J ^ 



where the coefficient C(s) is has been defined in equation G3.17|) . An important 
distinction of this equation is a new term with n = 0, that was absent in the equation 
(p.8|) for the Dirichlet boundary condition. 



We again use the uniform asymptotic expansion |29], ^] but now for the derivative 
of the Bessel function l v {yy) 

1 (1 + 1/ 2 ) 1/4 



KM 



V"/ 



V2 



y 



Vk(t) 



(3.28) 



k=l 



where rj is defined in equation (|3.10| ) and the functions Vk{t) are the known polynomials 
| 28fl in t(y). Again we use the approximation 

Vk(t) 



In l+£ 



k=l 



V 



E 

k=l 



F k N (t) 



where 



F 2 N (t) 
F 3 N (t) 
F A N (t) 



3 

1 



24 



--t 2 + -t 4 
16 8 

369 



16 



21 
128 



315 



1463 



r + —t° 1 h r, 



640 



128 32 



128 
733 

"64 



1152 
441 



32 



t 10 - 



707 
128 



12 



Instead of equation (|3.13j ) we have now 

4 



(3.29) 



(3.30) 



(3.31) 



where 




Z N (s) 


= C(s) [ 
Jo 




= Z^(s), 


Zf{s) 


= C(s)p- 



dyy l - 2s —HyI' (y)\ 



dy 



Z^s) , 



2s-l+j 



poo J 

{ R (2s+j-l)J o dyy^-F^t), j = 1,2,3,4. 



(3.32) 
(3.33) 
(3.34) 



1, . . . , 4 in the expansion ( |3.31D presents no 
difficulty while the function Z N (s) needs more thorough treatment. By making use of 



Finding the functions Z^(s), j 
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the asymptotics 

^-ln[yl' (y)]~ l + ly + 0(y 3 ), y^O, (3.35) 

^ yHylM] ^ 1 + ± + ± 2+ ^ + J^ + 0(y-% S ^oc (3.36) 

it is easy to make sure that the integral in the definition of the function Z (s) does not 
converge at any values of s. In order to overcome this drawback first the zeta function 
Cn{s) for a scalar field with nonzero mass m should be constructed and on calculating 
the residua according to equation (|3.1|) with this zeta function the mass m should be 
put equal to zerc0[|. Following this line we consider the function 

dy(y 2 -m 2 )-* +1 / 2 -HyI' (y)}. (3.37) 

It is defined in the region 1 < Re s < 3/2, the lower (upper) limit in these inequalities 
being determined by the convergence of the integral when y — > oo (y — > m). 

In order to find, by equation (^.1|), the contribution of the function Z^(s) to the 
coefficient B the analytical continuation of Z^(s) to the region Re s > 3/2 is required. 
In the most simple way it can be done by adding and subtracting from the integrand 
its asymptotics when y — > m. For our goals it is sufficient to take only the first term in 
this asymptotics 

= f°°r\ii (n 2 — m 2 ') 1 / 2 ~ s f — ]r,\iiT'JiA] — f \+f r\n In 2 — m 2 W 2 " 

COO 

where 



J°°dy (y 2 - m 2 ) 1 ' 2 -* |^ Hyl' (y)} - / m J + f m J°°dy (y 2 - m 2 ) 1 / 2 -, (3.38) 



fm= d\ ln[yI '° 



y=m 

The first term in equation (|3.38|) is regular at the point s = 3/2, but the integral in the 
second term gives a simple pole at this point 

ay[y ] 2 r(i/2) 

However the gamma function T(3/2 — s) responsible for this pole is canceled by the 
same multiplier in the denominator of the coefficient C(s) (see equation (|3.17j )). As a 
result the function Z^(s) does not give the contribution to the heat kernel coefficient 
B . 

For calculating the contribution of the function Z^{s) to the rest heat kernel 
coefficients by equation ( [3.1|) , this function should be analytically continued to the 
region Re s < 1. It can again be done by adding and subtracting the asymptotics 
of the integrand when y — > oo now. Further the residua are found according to 
equation ( |3.1| ) and only after that the mass m is put equal to zero. The corresponding 
results are presented in table 2. Finding the contributions of the rest functions 
Z^(s), = j — — 1, . . . , 4 to the heat kernel coefficients presents no trouble (see table 
2). 

Another way to treat n = case for the Neumann boundary conditions has been explained in detail 
in |3l. 
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3.3. Spectral problems on a cone 

Imposing on the eigenfunctions the periodicity condition with respect to the angular 
variable 9 with a period a we pass from a wedge to a cone because in this case the 
respective points on the radial rays 01 and 02 are identified (see figure 1), the circular 
arc 1-2 being converted into a circumference. Apparently, the boundary discontinuities 
at the points 1 and 2 disappear. This circumference separates two parts (internal and 
external) of the cone surface which for simplicity will be refereed to as the internal 
and external sectors of the cone. On the circumference separating them we impose on 
the eigenfunctions, as before, the Dirichlet or Neumann conditions. The corresponding 
(unnormalized) eigenfunctions for internal cone sector are 

u nm (r, 6) = Jnp(X nm r) ( smnp ® ) p = — n = 0, 1, 2, (3.39) 

\cosnp9 J a 

where X nm are the roots of the frequency equations for the Dirichlet boundary conditions 
at r = ill 

Jn P {KmR) = 0, n = 0, 1, 2, . . . (3.40) 

or of those for the Neumann conditions 

J' np (\ nm R) = 0, n = 0, 1, 2, . . . . (3.41) 

For the external sector the Bessel functions in equations (|3.39|) - (|3.41|) are replaced 
by the Hankel functions Hn P \\r). 

From ( |3.39|) it follows that on the cone all the states with n ^ are double 



degenerate 

N = l, N n = 2, n = 1, 2, . . . . (3.42) 

Thus, in order to find the heat kernel coefficients for the internal cone sector one 
should put in previous calculations 

p = 2n/a, (3.43) 

take into account the degeneracy of states (|3.43|) , and sum up over n starting with n — 0. 

Let us consider the internal sector of a cone with the Dirichlet conditions on the 
circumference 1-2. For the corresponding zeta function Cd( s ) ^he representation ( |3.8| ) 
holds with p defined in ( |3.43|) and with the summation replaced by 2 ^'^1 , where the 
prime on the summation sign means that the n = term is counted with half weight. 



For Cd( s ) m sum ( 3.13|) there arises an additional term with n = 



where 



C D (s)=Z D (s) + 2j2z?(s), (3.44) 

i=-i 

Z D (s) = C(s) / dy y 1 ^- In I (y) , (3.45) 

Jo d y 



and the functions Z^(s), j = —1, 0, ... ,4 are defined by equations (|3.14|) - (|3.16| ) with 
p = 2n/a. 
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The asymptotics 

±lnI (y)^^-^ + O(y 5 ), y^O, (3.46) 

-\nI (y)^l--- — - — --^ 1 + O(y- 5 ), y ^ oo (3.47) 
ay 2y 8y 2 8jr 12% 4 

imply that the function Z D (s) in equation ( 3.45Q is defined in the region 



1 < Re s < 3/2. (3.48) 
To single out in the integral (|3.45|) the pole contribution at the point s = 3/2 we rewrite 



the function Z D (s) as follows 

rDf\ nf„\ I J l-2s 



Z D {s) = C{s) f'dyy 1 
Jo 

+ dyy2 ' 2s + c(s) f dyyl ~ 2s a\ Io{y) ■ (3 - 49) 

The integrals in the first and third terms in equation ( |3.49| ) are regular at the point 
s — 3/2. Substituting the second term from equation ( |3.49j ) into definition ( ^.1| ) we get 

Res (Z D (s)T(s))= Res ( _L_")= . 

s^3/2-0 5^3/2-0^2^^(3/2 — 5)3 — 25/ 

Thus the function Z D (s) does not give any contribution to the coefficient B . 

So as to find the contribution of the function Z D (s) to the heat kernel coefficients 



B n , n = 1/2, 1, ... we again split the domain of integration in equation (|3.45| ) into 
two intervals (0, 1) and (1, oo). When integrating over the second interval we add 
and subtract under the integral sign the asymptotics ( |3.47| ). When calculating the 
residua at the points s — 1, s = 1/2, 0, —1/2, —1 we shall take the right-hand limits 
in equation ( p.l|) . It gives the following contributions to the heat kernel coefficients 
£1/2, B lt B 3/2 , B 2 , and B 5/2 , respectively 



2RyfR t -7T, -4. "^4- (3-50) 

In order to evaluate the contributions of the functions Z®(s), j = —1, 0, . . . , 4 to 
the heat kernel coefficients, in addition to (|3.50|) , one should substitute in table 1 a by 



a/2 and multiply all the elements of this table by 2. Summing the contributions (|3.50| ) 
and the data from table 1 we obtain the heat kernel coefficients for the internal sector 
I on the cone with the Dirichlet conditions on the circle 1-2 

aR 2 „ 1 „ ^ „ 2 n 2 a 



B = —-, B 1/2 = --aR^, Bi = - — + 
2 2 3 a 6 

a \/n „ 4 a „ 37 a \ fn 



B ' i/2 "64i?' ^"SIS^' B V 2 ~—2^- (3 ' 51) 
For Neumann boundary conditions the spectral zeta function for internal sector on 
a cone is given by 

4 

CM*) = 2 £zf( a ), (3.52) 
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where the functions Z N and are determined in equations ( |3.32| ) - ( |3.34| ) with 



a 



replaced by a/2. The corresponding coefficients of the heat kernel expansion are 

r, aR2 r, 1 j~y t — 2 7r 2 a 

^o = ^—, Bi/2 = - a R Vtt, B l = - — + - ) 
2 2 3 a 6 



_ 5a^ 4a _ 269 ay/? 

53/2 -"641"' B2 "45lP' ^-^I^- (3 ' 53) 



4. Identification of the individual contributions to the heat kernel 
coefficients for internal sector 

Let us envisage the internal sector I (see figure 1). Its boundary possesses the following 
peculiarities which give contribution to the heat kernel coefficients: nonzero curvature 
of the arc 1-2; right-angled corners at the points 1 and 2; corner of angle a at the 
origin. The arc contribution is proportional to its length, i.e., to a, the contributions of 
the right-angled corners does not depend on a, contribution of the corner at the origin 
vanishes when a = ir. It is sufficient to separate in each of the heat kernel coefficients 
the contribution due to each boundary singularity enumerated above. We demonstrate 
this considering the heat kernel coefficients for internal circular sector I with Dirichlet 
conditions on its boundary (see table 1). 
The first coefficient 

B = a^ = \n\ (4.1) 

is the area |0| of the circular sector I. The second heat kernel coefficient 

^2R+aR j— L 
B 1/2 = -Vn = -v^2' ( 42 ) 

where L is the length of the sector boundary (its perimeter). As concerns the coefficient 

the situation is more complicated. It is clear that ir/2 is the contribution of two right- 
angled corners at the points 1 and 2. The term a/6 contains the contribution of the 
curvature of the arc 1-2 (denote it by fc arc a) and a part of the contribution of the corner 
at the origin O (the latter is equal to (1/6 — /c arc ) a). In terms of these notations the 
complete contribution of the corner at the origin is 

c(a) = ^-+(\- k. Ar }\ a . (4.4) 
o a \6 / 

From the condition 

c(tt) = 7T Q - K T )j = (4.5) 

it follows that 

fcarc = 1/3 . (4.6) 
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Further we find c(a) 



n 2 -a 2 



c(a) = — . (4.7) 

o a 

Thus the corner of an angle a on the boundary gives the contribution to Bi defined 
by equation (|4.7| ). At first time this contribution has been calculated in [16]. Another 



method to derive it is described in [fL7P . In both the cases the Green's function of the 
equation of heat conductivity was considered. In our approach it is found by making 
use of the spectral zeta function technique. 

Now we are in position to check the consistency of our reasoning, namely, we can 
calculate the contribution to the coefficient Bi due to the right-angled corners at the 
points 1 and 2 by making use of equation Q4.71) . It gives 

2 c(tt/2) = vr/2 . (4.8) 

It is this value that has been attributed to this contribution above (see equation (|4.3|)). 

Identification of individual contributions to the rest of heat kernel coefficients can 
be done in a direct way. The terms independent of the angle a are attributed to the 
right-angled corners at the points 1 and 2, while the linear in a terms are due to the 
curvature of the arc 1-2 (see table 3). Merely such a separation of individual peculiarity 
contributions leads to a correct value of the arc curvature contribution which is known 
from the heat kernel expansion for a smooth boundary (see below). It is worth noting 
that the angle a at the origin does not contribute to the heat kernel coefficients B n with 
n > 1 even when a = tt/2. It may be explained only taking into account that the higher 
derivatives of the radius vector of the boundary curve behave in a different way at the 
origin and at the points 1 and 2. 

Let Q be a simply connected region of a plane bounded by a smooth curve T. 
For the heat kernel of the minus Laplace operator with the Dirichlet conditions on F 
the following asymptotic expansion holds when t — > (see, for example, references 

1,0,0]) 



37 / fc4(s)ds -2^/ (A;/(s))2ds 



+ 0(t 2 ) , (4.9) 



_2 15 7T 

where \Q\ is the area of Q, L is the length of T, k(s) is the curvature of the curve T at the 
point s, k 2 (s) = (d 2 r/ds 2 ) 2 , where r(s) is a parametric representation of the curve T; s 
is the natural parameter along T: ds 2 = (dr) 2 ; k'(s) = dk(s)/ds. For convex portions 
of T k(s) is considered to be positive, and for concave parts of T k(s) is assumed to be 
negative. 

In the expansion ( f4.9| ) the numerical coefficients of k, k 2 , k 3 , and k 4 are derived 
from the contributions, proportional to a, to the heat kernel coefficients Bi, -B3/2, B 2 , 
and -B5/2, respectively (see table 3). Here it should be taken into account that in the 
problem under study k(s) = l/R, ds = Rda and the coefficients B n enter the heat 
kernel expansion (|2.5| ) with the multiplier 1/(4 7r). 
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If we go from a wedge to a cone by identifying the radial rays 01 and 02 the 
corners at the points 1 and 2 disappear. The sole singular point remains the origin O 
which becomes an internal point of the cone surface. However the contribution of this 
singularity to the coefficient Bi has the same nature as in the case of wedge. 

The heat kernel coefficients for internal sector I on the cone with Dirichlet condition 
on the boundary 1-2 are listen in equation (|3.51 ). The coefficients B and By 2 obey 
general equations (|4.1|) and ( [4.2|) with corresponding values for and L. Separation 
of the contributions to B\ generated by singularity at the origin and by the curvature 
of the circle 1-2 can be conducted in the same way as for a wedge. Let k arc a be the 
contribution of the boundary 1-2 and 

2 7r 2 a 

d(a) = - — + --h Ic a (4.10) 

3 a o 

be the contribution due to the singularity at the origin. When a = 2ir the surface 
of a cone becomes a plane and the singularity at the origin O disappear. Therefore, 
d(2n) = 0. It gives 

fcarc = l/3, (4.11) 

i.e., for this quantity we have the same value as in the case of a wedge (see (|4.6|) ). With 
Ql.llp allowed for, one deduces from the equation ( |4.10|) 

d(a) = |— -^ = 2c(a/2), (4.12) 
S a o 

where c(a) is defined in equation (|4.7|) . 

Identification of the individual contributions to the heat kernel coefficients for the 
Neumann boundary conditions is conducted in the same way (see table 3). The first 
coefficient B is, as before, the area \Q\ of the sector I. For the second coefficient By 2 
wehave equation fl4.2|) with oposite sign in the right-hand side 

L 
2' 

The contribution to the coefficient B\ generated by the singularity at the origin proves 
to be the same as for the Dirichlet conditions, c(a). 

For Neumann boundary conditions we failed to find in the literature the formula 
analogous to equation (|4.9| ), i.e., the asymptotic expansion of the heat kernel for the 
operator —A defined in the region of a plane with smooth boundary curve. However, as 
was noted earlier, the spectral problem on the plane, envisaged by us, and its cylindrical 
generalization have the same coefficients B n . Therefore, for verification of our results, 
concerning the contributions to B n due to the smooth parts of the boundary, we used 
the expansion of the heat kernel for Robin conditions with smooth boundary that has 
the dimension greater than 1 (see references [0, 0)- 



B 1/2 = v ^- } L = 2 R + a R. (4.13) 
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5.1. Basic formulas 



Here we shall use the technique applied in ||23|| . It is close to the method developed in 
@, 2(| |2l] and employed in preceding subsections. 

Let us consider the spectral zeta function depending on a parameter x 2 



It may be regarded as an extension to the general spectral problem of the Epstein- 
Hurwitz zeta function 



CEH<>,a 2 ) = + 



n=l 



It turns out that the heat kernel coefficients B n can be found from the expansion 
of the function £(s, x 2 ) in terms of inverse powers of x developed for a certain value of 
s. It is convenient to chose this value to be equal to d/2. In fact, from the definition of 
the gamma function it follows that 



dtt s ~ l e- {x2+x2)t . 



T(s) (X 2 + x 2 )- s = 
For s = 1 + d/2 equation ( |5.2|) gives 



(5.2) 



-xlt 



POO 

/ dtt d/2 e~ x2t K(t). (5.3) 
Jo 



On substituting the asymptotic expansion (|2.5|) in equation ( |5.3j ) we obtain 



r(l + rf/2)C(l + rf/2,x 2 )^ J2 



B 



n/2 



n=0 



( 4vr )rf/2 



r i + 



n 



x 



-n-2 



B | £?i/ 2 r(3/2) | B 1 T(2) | ^ 3/2 r(5/2) 



x- 



X' 



X' 



X" 



, ^r(3) + B 5/2 T{7/2) + B 3 T(4) + 



( 5 - 4 ) 

tAJ tXj %KJ 

Let F(z) = be the frequency equation which determines the spectrum A n in the 
problem under consideration. We also suppose that the function F(z) allows one to 
rewrite this equation in the form 

n( A n-^)=°- ( 5 - 5 ) 



Taking into account the equality 

1 (-l) m 



(Xl + x 2 ) m r(m) \2xdx 
we recast the left-hand side of equation 



ln(A 2 + x 1 



to the form 



z = ix . 



d 



d 



, 1 d\ l+d/ \ . 

( — \nF(ix). 

1 2xdxJ v ' 



(5.6) 



(5.7) 
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Obviously formula ( |5.7p is applicable only to the manifolds of even dimension. 

Rather than to calculate the heat kernel coefficients B^ 1 for the external sector 
alone it is simpler to find first the coefficients B^ +n for the union of the sectors I and 
II. Then the coefficients B^ 1 are obtained as the corresponding differences 

B» = Bi +n -Bl n = 0, 1/2, 1, . . . . 

When calculating the heat kernel coefficients for the union of the internal (I) and 
external (II) sectors on should take into account the following. The technique used 



by us U, [20|, |21], |23], |24[ gives the difference between the zeta function for the region 
I+II and the zeta function for the corresponding part of the Euclidean space. The last 
contribution is usually referred to as the Minkowski space-time contribution f33fl . In 
the case under consideration the zeta function for an open angle a (without circular 
arc 1-2) is subtracted from the zeta function sought for. The heat kernel coefficients 
B n , corresponding to the Minkowski space-time contribution, can be calculated in the 
following way: in the respective heat kernel coefficients for the sector I one should 
put R = R\ —>■ oo and discard the curvature contribution of the arc 1-2 and of two 
rightangled internal corners at the points 1 and 2. 

5.2. Internal and external circular sectors of a wedge with Dirichlet condition on 
separating arc 



Now we proceed to practical using the general formulas (|5.4|) and (|5.7p for calculation 



of the heat kernel coefficients. We consider the scalar Laplace operator on the union of 
internal and external circular sectors on a plane (see figure 1) with Dirichlet conditions 
on the arc 1-2 separating these sectors. As it was explained in Section 2 the heat 
kernel coefficients for the corresponding boundary value problem in space (d = 3) The 
frequency equations are will be the same. 

J np (Rz) = 0, (internal sector I), (5.8) 
H$(Rz) = 0, (external sector II), n = 1, 2, . . . , p = n/a. (5.9) 

Further we shall concern with the product J u (ix) He ,{ix) and use for it the uniform 
asymptotic expansion |2£| which depends only on x 2 . In this case the condition ( |5.5|) is 
apparently satisfied. This point can also be explained in the following way. In view of 
the formula |30, |34 



for the function z~ v J v (z) the representation of type ( |5.5| ) holds. Here z v ^ m (m = 
1, 2, 3, . . .) are the nonzero roots of the function J u (z). The multiplier {z/2) v in 
equation ( |5.10| ) is canceled in the product J u (ix)H^\ix) ~ I u (x)K v (x) with the 
multiplier (x/2)~ u following from the small x asymptotics of the function K v {x): 
K v {x) ~ (l/2)r(z/)(x/2) v , v > 0. Hence the requirement (|5. 5|) is satisfied. 
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Setting in fl5.7|) d = 2 and substituting in it the left-hand sides of the frequency 
equations (|5.8j) and (|5.9|) we obtain 

2 co 



(1 d \ 
' n=l 



(5.11) 



functions [28 



Now we use the uniform asymptotic expansion for the product of the modified Bessel 
\nI v (Rx)K v (Rx) -- 



ln2z/ + lnt + ^ 2jW 



2.; 



where 



1 



TIT 



v 



v = np 



n- 



a 



(5.12) 



(5.13) 



and G®{t) are the polynomials in t expressed in terms of the known functions Uk(t). In 
order to calculate the first six coefficients in the expansion (|5.4j ) it is sufficient to keep 
two terms of the sum in (|5.12|) . The relevant coefficients G%{t) and Gf(t) are 



4 



13 



-r 



71 



531 



339 



t 10 + 



565 



-t 



12 



(5.14) 
(5.15) 



64 16 32 16 64 
Substituting in (|5.11|) differentiation with respect to x by differentiation with respect to 
t we get 



Cd(2,x 2 



n=l 



R 



i 3 d 



n=l 



2 dt 

2^4 + 4^^ 



lnt + ^Gg(t) 
j'=i 



18t 



=i 

1 ' 30t 4 ) 



3t 8 



32z/ 8 



;i3 - 568t 2 + 3540t 4 - 6780t 6 + 3955t 8 ) 



(5.16) 



It will be recalled that t depends on x and n through ( 5.13 ). 

All the sums in the equation (|5.16 ) are finite. Hence the problem of analytic 
continuation does not emerge here. In order to do the summation exactly we use the 



formula 30, 35 



E 

n=l 



rr 



7T 

2y V 



( coth 



ny 



ny 



7T 



1 



2y 2y 2 



Si(y). 



(5.17) 



When y — > oo, the function affords the value of the sum on the left-hand side of 

equation ( |5.17| ) up to exponentially small corrections. Step-by-step differentiation with 
respect of y of the left and right hand sides of the equation ( |5.17|) gives the values of all 
the sums entering (|5.16|) 
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00 1 X X d 37r X 



1 1 d 



Sh(y) = S k +i(y) 



(5.18) 



~[ (y 2 + n 2 ) k+1 k 2y dy 

In order to express the zeta function ( D (2,x 2 ) in terms of Sk(x), k = 1, 2, ... ,8 
explicitly the following substitutions should be accomplished in equation fl5.16| ) 

— - 1 

,.2 ~~ -9. 1 ^?.\ ' (O.lyj 



p 2 [n 2 + x 



x 2 



n 2 + x 2 



x 



Rx 
P 



(5.20) 



the change ( 5.2(J| ) being done only in the round brackets in this equation. As a result 
the zeta function ( D (2,x 2 ) acquires the form 



Cd(2,x 2 



$ S2 ^)-$[^ S ^) 



42x 2 S 4 (x) + 30x i S 5 (x 



3R 4 



[160 S. 



AX 



1992x 2 S 5 (x) + 6930x 4 S 6 (x 



32p 8 

-9040 x 6 S 7 (x) + 3955 x s S s (x)] . (5.21) 

Substitution of the explicit expressions for the functions Sk(x) from ( |5. 18 ) to ( |5.21| ) 
gives 



Cd(2,x 2 )= -^ + £t4 



3 a 



+ 



+ 



555 a 



+ 



39 



ix° 



4x 4 512i?x 5 8R 2 x 6 2 17 R 3 x 7 2 6 i? 4 x 8 



+ 0{x 9 



(5.22) 



Comparing the expansions (|5.4j ) and (|5.22 ) we obtain the values of the first seven heat 
kernel coefficients in the problem under consideration 



BI+II 



B 



i+n 



B I+ 

-°3/2 



0, ±jy 2 

11 ay/n 



-aRy/n, B[ 



i+ii 



n 2 



7T 



51+11 
n 5/2 



7T, 

31 \pn a 



51+11 
n 3 



13tt 



(5.23) 



32 R' " z 4i? 2 ' 5 / 2 2 12 R 3 ' ~ 6 32R* 
where B^ +n = B^ +n — B n , and B n are the heat kernel coefficients corresponding to the 
open angle a (without the arc 1-2) 



B = - aR 2 , B 1/2 = -y/wRx 



aRi, By 



ir 2 -a 2 



clot) 



Ri 



00 . 



(5.24) 



All the rest of coefficients Bj vanish Bj 



6a 

0, j>3/2. 

Using equations ([5.23 ), ( 5.24 ) and the results for the internal sector I, presented in 
table 1 we deduce the heat kernel coefficients for the external sector II with Dirichlet 
boundary condition 



B 



11 



R 2 



711 



B l/2 



-V^(Ri-R) 



B 



11 

3/2 



7r a \/7r 
AR + 64 i? ' 



B. 



11 



71 



4 







8R 2 315 R 2 



^(R 1 + R), B[ T = — - — 
2 v it j, 1 2 3 

25 -Jit 37 a t/tt 
~96lV + 2 13 R 3 ' 



^ B 



11 

5/2 



(5.25) 
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5.3. Internal and external circular sectors of a wedge with Neumann condition on 
separating arc 

In the case of Neumann boundary conditions on a separating arc the heat kernel 
coefficients for the union of internal and external sectors are calculated completely in a 
similar line. Now the frequency equations read 



0, (internal sector I), (5.26) 



rz 



T=R 



r=R 



0, (external sector II), n — 0, 1, 2, . . . , p = rr/a. (5.27) 



The angular part of the eigenfunctions is proportional to cos(np6), therefore the index 
n takes integer values starting with zero. 

For the zeta function in this eigenvalue problem the representation analogous to 
( CT ) holds 

/ 1 A \ 2 00 

Cn(2,x 2 ) 

71=0 

Again we use the uniform asymptotic expansion 



1 d \ 2 °° 

/ n— n 



(5.28) 



In 



ln2^-lnt + ]T^ (t) 



Ai 



(5.29) 



3=1 



where 



G N(A = --t 2 + -t 4 - -t\ 

2 w 8 4 8 ' 



G l l{t) 



109 



27 4 

— t 4 + 

64 16 



441 



733 „ 

1° + — 

32 16 



t 10 - 



707 

i 

64 



12 



(5.30) 



On substituting (|5.29|) and (|5.30|) and differentiating in equation ( |5.28 ), the zeta function 
under consideration assumes the form 



0v(2,x ! 



n=0 



t 4 



3t e 



2z/ 4 4z/ 6 



'1 - 10t 2 + 14 1 2 ) 



32zv 8 



'81 - 2616 1 2 + 14660 1 4 - 26460 1 6 + 14847 1 8 ) 



(5.31) 



In order to take the sum over n in (|5.31 ) we apply the formula that follows from ( 5.17|) 



oo 1 
x 1 7T 

' j/ 2 + n 2 2y 



n=0 



ny 



71 



coth7ry H ~ — + — = Si(y) . 



2y 2^/2 



Sequential differentiation of equation ( |5.32|) gives 

1 lid 



E 

n=0 



(j/ 2 + n 



k 2y dy 



Sk(y) = S k+X (y) 



(5.32) 



(5.33) 
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On making use of the change of variables ( |5.19| ) and ( |5.20| ) in equation ( |5.31| ) the zeta 
function (jv(2,x 2 ) assumes the form 

(n(2,x 2 ) = ^S 2 (x) + ^ [ 5 &(*) ~ 18x 2 5 4 (x) + Ux 4 S 5 (x)} 

+ [512 ^(x) -6712 x 2 S 5 (x) + 24362 x 4 S 6 (x) 

32p s L 

-32928 x 6 S 7 (x) + 14847 x s S 8 (x)] . (5.34) 

Substitution of the function Sk{x) with k = 2, 3, ... ,8 from (|5.32j) and (|5.33|) gives 

, , 2\ ql R 1 15 a 3 1 4035 a 81 1 . _ 9 . , . 

CM ,X j = 8^ + 4^ + 512 Wx^ + 8 ^ + 131072 RW + 64 Rh^ + {X j ' ( } 
Comparison of the expansions ( p.35| ) and ( |5.4| ) gives the following values for the heat 
kernel coefficients B^ +n for Neumann boundary conditions 

B 1 ^ 11 = 0, 5*+" = aRy^, B[ +n = vr, 

3 / 2 32 R ' 2 4# 2 ' 5 / 2 4096 R 3 ' 3 32 i? 4 y J 

The heat kernel coefficients corresponding to the Minkowski space-time contribution 
in the case under consideration are derived from the respective coefficients for internal 
sector I by putting there R = R± — > oo and omitting the contribution due to the 
curvature of the arc 1-2 and contributions of the rightangled corners at the points 1 and 
2 (see Tables 2 and 3) 

B = ^R 2 , B 1/2 = ^R 1 + ^^R 1 B 1 = c{a), B J = 0, j = 3/2, 2, . . . . (5.37) 



By making use of the equations Q5.36j ), (|5.37|) and table 2 we derive the heat kernel 



coefficients for the external sector II with Neumann condition 

Bi' = \ (Ri-R 2 ), B(j 2 = V^(Ri-R) + ^(Ri + R) *" = £-f 



u 3Vtt | 5a0r u _ 3tt 4 a n _ 21 A 269 a yfn 
^ 2 ~ AR + 64i? ' * 2 ~8R 2 45 R 2 ' 5/2 " 32 R* + 2^3 'V b - d8 J 

5.^. Internal and external sectors on a cone 

In the case of a cone the eigenfunctions are defined in ( |3.39| ) with X nm being the roots 
of the frequency equations (|5.8|), ( |5.9j ) and ( [5.26D , ( |5.27| ) where p = 2n/a. For both 



Dirichlet and Neumann conditions on the circle 1-2 the index n ranges from n = 0. The 
state degeneracy in both the cases is determined by equation ( |3.42| ). 



All this implies that in order to proceed to a cone one should put in the relevant 
formulas for a wedge p = 2-nja and use for summation over n the following relations 
(instead of equations (|5.18|) , (|5.32|) , and (|5.33|) ) 



1 oo 



— cotli7ry ~ — = S{(y) 



y + n y y 

^^t^^^-l—^y^ 3 ^- * = LV... (5.39) 
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It is essential that these summation formulas should be employed for both Dirichlet and 
Neumann boundary conditions given on the circle 1-2 and separating the internal and 
external sectors of a cone. 

First we consider Dirichlet conditions. Carrying out in (|5.21 ) the change of variables 
( |3.43| ) and using for summation the functions S%(x) from equation ( |5.39| ) we obtain the 
expansion for the zeta function in the spectral problem at hand 

.„ . 9 . aR 3 a 555 a Q , , > 

^■^-s?^ + WJ^ + °( X "')■ < 5 ' 4 °) 

This expansion can formally be derived from the relevant zeta function for a wedge, 
equation (|5.22|) , by omitting the terms with even powers of x. Comparison of the series 
( |5.40| ) with equation ( |5.4j ) affords the values of the first seven heat kernel coefficients 

n 

B^ 11 = B[ +n = Bi +H = Bi +n = 0, Bi+" = -a V^R, 
5i+n _ bi+ii _ 37a^ 

B ^ ~ 32R ' ~ 2«J2> • [bAi) 

The heat kernel coefficients for Minkowski spacetime contribution are 

B =-Rl B 1/2 = --^R 1 , B 1 = 2c(a/2), B, = 0, j > 3/2 . (5.42) 

We deduce from equations ( |5.41| ), ( |5.42j ), and ( p.51| ) the heat kernel coefficients for 
external sector of a cone with Dirichlet boundary conditions on the circle 1-2 

b q = 2 ~ b i/2 = -^Vn (Ri- R), B[ t = --, 



n _ a y/Tr jj 4a n _37ay/w 

64iT ^ ~ 315 i? 2 ' 5 / 2 " 2«ff ' V A6 > 



By analogy with this we obtain from equation (|5.34 ) for Neumann boundary 
conditions 

9 . aR 15a 4035 a q . . 

ftft , 2 ) = — 3 + + -p- + 0( X ■) . (5.44) 

This expansion can formally be derived from the relevant zeta function for a wedge 
equation ( |5.35| ) by omitting the terms with even powers of x. Comparing the equations 
( |5.44j ) and ( J5.4j ) we obtain the values of the first seven heat kernel coefficients B^ +n 

Bp" = B{ +n = B^ 11 = Bl +n = 0, Bi+ n = a yfrR, 

5/+I7 _ 5ay^ = /+// _ 269 a 

^ 3 / 2 " 32R ' ^ 5 / 2 " 2W ' ^-45) 
The heat kernel coefficients for Minkowski spacetime contribution in this case are 

B = ^Rj, S 1/2 = |0Fi2i, Bi = 2c(a/2), ^ = 0, j > 3/2 . (5.46) 

From equations ( 5.45|) , (|5.46 ), and ( |3.53|) we derive the heat kernel coefficients for 
the external sector II of the cone with Neumann boundary conditions on the circle 1-2 



(// 5a0r jj 4 a TI 269a^/7^ 



B3 / 2 ~ ~64iT' B2 "45JP' ^Z 2 " 2^3 • ( 5 - 47 ) 
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6. Specification of individual contributions to heat kernel coefficients for 
external sector. Some general rules 



The heat kernel coefficients for external sector are given in equations ( 5.25 ), ( |5.38| 



( |5.43| ), and ( |5.47| ). Identification of the individual contributions to these coefficients due 
to the boundary nonsmoothness does not differ basically from the analogous procedure 
for internal sector carried out in section 4. The sole point that should be noted here is 
the following. The area and the perimeter L are in this case infinite. Therefore they 
should be treated as the limit, when Ri tends to infinity, of the expressions 

|0| = ^{R\-R 2 ), L = 2 — R) + a(Ri + R) . (6.1) 

The coefficients B and Bi/ 2 are defined, as before, by equations (fP]) , ft4.2|) , and (|4.13| ) 
with \Q\ and L defined in ( |6.1|) . Contributions to the rest of the heat kerenel coefficients, 
which are proportional to the angle a are due to the curvature of the arc 1-2 and the 
contributions which are independent of a are due to the rght-angled corners at the 
points 1 and 2. The results of this analysis are represented in table 3 together with the 
heat kernel coefficients for the internal sector I. Analysis of this table enables one to 
reveal some general rules for corner contributions. 

It is interesting to note that the right-angled corners at the points 1 and 2 give the 
contributions to the coefficients B 3 / 2 and B 5 / 2 which have opposite signs for internal and 
external sectors, it being valid for all configurations and boundary conditions considered. 
Such a behaviour seems to be related with convexity (internal sector) or concavity 
(external sector) of the arc 1-2. It is analogous to the contributions of the smooth 
segments of the boundary to the heat kernel coefficients B\ and B 2 . 

The boundary discontinuities due to the corner at the origin O and to the corners 
at the points 1 and 2 contribute to the heat kernel coefficients basically in different ways. 
Corner at the origin gives contribution only to the coefficient Bi (even when a = 7r/2). 
The corners at the points 1 and 2 contribute to all the heat kernel coefficients starting 
with n = 1. Obviously the reason of this distinction is a geometrical one, the corner at 
the origin is formed by crossing two straight lines, while the corners at the points 1 and 
2 are the result of intersection of lines one of which has a nonzero curvature. 

This general assertion concerning the corner contribution to the heat kernel 
coefficients can be illustrated by a known expression for the heat kernel expansion for a 
rectangle with sides a and b 



OO OO , ., ., N , 1 \ / 7 1 

2 . I a 1\ / b 1 



m=l n=l 



2 2 

m n 



Ant 2 J V JAnt 2 



ES 



Ant 4v^rt 4 

where ES denotes the exponentially small corrections as t — > (see, for example, ||19|| ). 
Here the scalar operator —A with Dirichlet boundary conditions is considered. For a 
rectangle the coefficient B\ is obviously equal to the contributions of four right-angled 
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corners. Indeed, the third term in the expansion (|6.2j) can be represented in the form 

11 1 

- = — B x = — 4 da = vr/2) , 
4 4vr 1 4vr V 1 



where c(a) is given in equation (|4.7|) . Besides these corners the boundary of a rectangle 
has no other singularities, therefore the heat kernel coefficients B n with n > 1 vanish in 
this problem. 

These rules for obtaining the heat kernel coefficients are directly generalized to an 
arbitrary polygon with the angles a.i. The first two coefficients Bq and B1/2 are defined 



by equations (4TT) and ( |4.2j ), respectively, where is the area of the polygon and L is 
its perimeter. The third coefficient Bi is equal to the sum of the contributions due to 
the angles «j 

B 1 = Y^<<x i ). (6.3) 

i 

The rest of the coefficients B n , n > 1 vanish. In particular, it implies that the zeta 
function technique should provide a finite value of the Casimir energy for a polygon on 
a plane (B 3 / 2 = 0) and for a cylindrical generalization of the polygon spectral problem 
(B 2 = 0). These subjects have been discussed earlier in papers [P^j . 

In reference |J7| the vacuum energy of massless fields including electromagnetic 
field was calculated for the boundary configuration shown in figure 1 with a = ix. Both 
versions of this boundary value problem were considered, three-dimensional one (a semi- 
circular infinite cylinder) and two-dimensional spectral problem on the plane. In both 
the cases the zeta regularization didn't give a finite value of the Casimir energy. As 
known ||, the reason of this is nonzero heat kernel coefficients B 2 for d = 3 and B 3 / 2 
for d = 2. Using the results of the present paper we can elucidate the geometrical origin 
of this fact. Let us consider electromagnetic filed in internal and external sectors of the 
wedge with a = it. For the boundaries with cylindrical symmetry electromagnetic 
field reduces to two massless fields subjected to Dirichlet and Neumann boundary 
conditions |38, 39| . From table 3 it follows that 

The nonzero value of the coefficient B 2 is due to the contribution of four right-angled 



corners at the points 1 and 2. It has been noted at first time in reference ||15|| . In the 
case of the coefficient B3/2 the contributions of the corners from internal and external 
sectors are mutually canceled, while the contributions of the curvature of the arc 1-2 
from internal and external sectors are added. 

Different geometrical origins of the zeta function failure to provide a finite value of 
the vacuum energy in the two- and three-dimensional versions of the boundary value 
problem in question probably imply the impossibility of obtaining a finite and unique 
value of this quantity by taking advantage of the atomic structure of the boundary EO] 
or its quantum fluctuations 4"1|. It is clear because any physical reason of the Casimir 
energy divergences should be valid simultaneously in the two- and three-dimensional 
versions of the boundary configuration under consideration. 
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7. Conclusion 

The basic result of the paper is the calculation of the individual contributions to the heat 
kernel coefficients generated by such particularities of the boundary as the corners. In 
the course of this analysis certain patterns, that are followed by these contributions, have 
been revealed. As a by product, the rules for obtaining all the heat kernel coefficients for 
the minus Laplace operator defined on a polygon or in its cylindrical generalization are 
formulated, these rules being valid both for Dirichlet and Neumann boundary conditions. 
Implications of the obtained results in calculations of the vacuum energy for regions with 
nonsmooth boundary are discussed. Our calculations comport the conventional point of 
view according to which the heat kernel coefficients are determined, in the case under 
consideration, by the local properties of the boundary. 

In any case, the heat kernel coefficients obtained in the present paper can be used 
for verification of the general methods of calculating the contributions of boundary 
discontinuities to the heat kernel coefficients which may be developed in the future. 
These general methods must in particular allow one to calculate the contribution of an 
arbitrary corner to the coefficients B n with n > 1 in terms of the jumps of the derivatives 
of the boundary curve (or its geometrical invariants) at this corner, i.e., the formulas 
analogous to equation ( f4.7| ) should be found for B n with n > 1. 
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Figure 1. The cross section of a dihedral angle with circular boundary of radius R 
inside. 



Table 1. The contributions of the functions Z®, j = — 1,0, . . .4 to the heat kernel 
coefficients (internal circular sector I with Dirichlct boundary conditions). 
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Table 2. The contributions of the functions Z N and Z^ , j = —1,0. . .4, to the heat 
kernel coefficients (internal circular sector I with Neumann boundary conditions). 
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Table 3. The contributions of different parts of the boundary to heat kernel 
coefficients; D and N stand for the Dirichlct and Neumann boundary conditions for 
a wedge, Dc and Nc denote these conditions for a cone; the upper (lower) sign is 
referred to the internal I (external II) sector. 
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